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Abstract 

We describe two new species of springsnails (genus Pyrgulopsis) for populations from the middle Fork and 
upper East Fork of the Gila River Basin (New Mexico) that had been previously identified as P. gilae. We 
also restrict P. gilae to its originally circumscribed geographic range which consists of a short reach of the 
East Fork Gila River and a single spring along the Gila River (below the East Fork confluence). These 
thiee species form genetically distinct lineages that differ from each other by 3.9-6.3% for mtCOI and 
3.7-8.7% for mtNDI (the latter data were newly obtained for this study), and are diagnosable by shell 
and penial characters. Collectively the three species form a strongly supported clade that is distinguished 
from other congeners by the unique presence of two glandular strips on the dorsal surface of the penial 
filament. These findings suggest that the conservation status oi P. gilae, which was recently removed from 
the list of candidates for listing as endangered or threatened by the United States Fish and Wildlife Service, 
should be revisited and that the two new species may also merit protective measures given their narrow 
geographic ranges. 
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Introduction 

Pyrgukpsis is a large genus (137 species; Hershler et al. 2013) of freshwater gastropods 
that is distributed in North America west of the Mississippi River basin. The tiny species 
in this genus Uve in spring-fed habitats and usually have very small geographic ranges. 
Pyrgukpsis is a current focus of conservation efforts owing to threats posed by groundwa- 
ter extraction, livestock grazing and other anthropogenic activities (Hershler et al. 2014). 
Recent molecular studies have shown that several congeners are composites of genetically 
divergent lineages and are in need of taxonomic revision (e.g., Liu et al. 2003, Liu and 
Hershler 2007, Hershler and Liu 2010). This is the second in an anticipated series of 
papers that clarifies the taxonomy of these species (Hershler et al. 2013). 

Pyrgulopsis gilae (Taylor, 1987) was described for specimens from single springs 
along the lower East Fork (type locality) and main stem Gila River in Grant County, 
New Mexico. Field surveys in the 1990's and 2000's resulted in the discovery of new 
populations in two other reaches of the upper Gila River watershed (Middle Fork, upper 
East Fork) that are currently being treated as P. gilae (NMDGF 2012). Hurt (2004) de- 
lineated substantial divergence in mtCOI sequences (6.8% average) between specimens 
from the upper East Fork reach (Wall Spring) and three localities within the originally 
circumscribed range of P. gilae. In a more comprehensive survey of COI variation within 
P. gilae, populations from the upper East Fork, Middle Fork, and lower East Fork (and 
main stem Gila River) reaches were resolved as three divergent (3.9-6.3%) sub-clades 
which were postulated to be distinct species (Liu et al. 2013). Here we document a con- 
gruent pattern of variation in a second mitochondrial DNA marker (NDI) and delineate 
morphological differences supporting recognition of the upper East Fork and Middle 
Fork Gila River populations as new species, which are described herein. 



Methods 

For the current molecular study we used the same samples that were analyzed in our pre- 
vious phylogeographic investigation of P gilae across its entire geographic range (Liu et 
al. 2013; Fig. 1). Genomic DNA was extracted from single, entire snails using a CTAB 
protocol (Bucklin 1992); 3-8 specimens were analyzed (separately) from each sample. 
ND43F and RND592F (Liu et al. 2003) were used to amplify a 530 bp fragment of 
NADH dehydrogenase subunit I (NDI). This primer pair did not amplify the region for 
specimens from two localities (G8, Gil) and consequently we designed a second set of 
oligonucleotide primers for these snails, ND30 (5'TCTTAYATR CAR ATW CGT AAA 
GG3') and RND490 (5 ATG TTA CAA ATC ATATAA ATG3'), based on conserved re- 
gions of NDI in an alignment from P gilae and two closely related species {P deserta [Pils- 
bry], P. davisi [Taylor]). Degenerate positions are represented by the following ambiguity 
codes: Y=C/T; R=A/G; W=A/T. Amplification conditions and sequencing of amplified 
polymerase chain reaction product followed Liu et al. (2003). Sequences were determined 
for both strands and then edited and aligned using SEQUENCHER™ version 5.0.1. 
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Figure I . Map showing the distribution of mtDNA clades I-III. 



The 56 newly sequenced specimens of 7? gilae were analyzed both separately and together 
with our previously published COI dataset (Liu et al. 2013). We included the same set of 
outgroup taxa as in our prior study of B gilae (Liu et al. 2013), with Fbridobia floridana 
again being used as the root. The GenBank accession numbers for these sequences (COI, 
NDI) are given in Appendix 1 . Note that we newly sequenced specimens of P. ''mimbres" 
for COI (using the methods of Liu et al. 2003) and NDI as part of this study (GenBank 
accession numbers: COI, KM205358; NDI, KM205359). The newly obtained haplo- 
types from ezchU gilae sampling locality were deposited in GenBank. Sample information 
and GenBank accession numbers are given in Table 1. One example of each haplotype 
detected in a given sample was used in our analyses. 
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Table I . Sample codes, collection localities, and GenBank accession numbers for /? giUe mtDNA sequences. 



Code 


Locality (all in New Mexico) 


COI 


NDI 


Gl 


Spring along East Fork Gila River, ca. 1.53 km north, 2.9 km east of 
State Route 527 bridge crossing, Grant County 


KC571284, 
KC571285 


KM079175, 
KM079176 


G2 


Spiing along East Fork Gila River, ca. 1.29 km north, 0.56 km west of 
Black Canyon confluence, Gtant County 


KC571286, 
KC571287 


KM079177, 
KM079178 


G3 


Spiing along East Fork Gila River, ca. 1.53 km north, 2.38 km east of 
State Route 527 bridge crossing, Grant County 


KC571288, 
KC571289, 
KC571290 


KM079176 


G4 


Seepage along Taylot Cteek, ca. 0.32 km south, 0.93 km west of Wall 
Lake dam (below Wall Lake), Catron County 


KC571291, 
KC571292, 
KC571293, 
KC571294 


KM079180 


G5 


Hillside seep, 1.61 km noith, 0.97 km east of Burnt Corral Canyon, 
Catron County 


KC571295 


KM079181 


G6 


Spring along Beaver Creek, ca. 0.29 km north, 0.40 km west of Taylor 
Creek confluence, Catron County 


KC571296, 
KC571297, 
KC571298 


KM079181, 
KM079182 


G7 


Seepage along Taylor Creek, 50 m west of Whitetail Canyon, Catron 
County 


KC571292 


KM079183 


G8 


Spiing along Middle Fork Gila River, ca. 0.97 km north, 0.64 km west 
of Jordan Canyon, Grant County 


KC571299, 
KC571300, 
KC571301, 
KC571302 


KM079187, 
KM079188 


G9 


Fall Spring, 1.61 km north, 0.56 km east of Burnt Coital Canyon, 
Catron County 


KC571303 
KC571304 


KM079184, 
KM079185 


tGlO 


Fall Spring, 1.61 km north, 0.56 km east of Burnt Coital Canyon, 
Catron County 


KC571304 


KM079185, 
KM079186 


Gil 


Spring along Middle Fork Gila River, ca. 0.48 km north, 0.48 km west 
of Jordan Canyon, Grant County 


KC571305 


KM079187, 
KM079189, 
KM079190 


G13 


Spiing alongTaylor Creek, 0.81 km north, 1.13 km east of Wall Lake 
Dam, Cation County 


KC571292 


KM079183 


G14 


"Alum Hot Spring," ca. 1.93 km south, 0.16 km west of State Route 
527 bridge crossing, Grant County 


KC571288, 
KC571306 


KM079179 



tVery small (juvenile) specimens initially thought to be distinct from P. gilae. 



The partition homogeneity/incongruence length difference test (F arris et al. 1994) 
was used to determine whether the COI and NDI datasets were consistent and could 
be combined for the phylogenetic analysis. This test, which was conducted using parsi- 
mony-informative sites only and 1,000 replicates, did not detect significant incongru- 
ence {P=0.2\) and consequently we combined the two datasets in the phylogenetic 
analysis. MRMODELTEST 2.3 (Nylander 2004) was used to obtain an appropriate 
substitution model (using the Akaike Information Criterion) and parameter values for 
this analysis. Phylogenetic relationships were inferred by Bayesian analysis using MR- 
BAYES 3.1.2 (Huelsenbeck and Ronquist 2001). Metropolis-coupled Markov chain 
Monte Carlo simulations were run with four chains (using the model selected through 
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MRJMODELTEST) for 2,000,000 generations, and Markov chains were sampled at 
intervals of 10 generations to obtain 200,000 sample points. We used the defauk set- 
tings for the priors on topologies and the GTR + I + G model parameters selected by 
MRMODELTEST as the best fit model. The Tracer program was used to analyze runs 
for Effective Sample Size (ESS, greater than 200) to ensure that sufficient sampling 
occurred. At the end of the analysis, the average standard deviation of split frequen- 
cies was less than 0.01 (0.002) and the Potential Scale Reduction Factor (PSRF) was 
1, indicating that the runs had reached convergence. The sampled trees with branch 
lengths were used to generate a 50% majority rule consensus tree with the first 25% of 
the samples removed to ensure that the chain sampled a stationary portion. 

Genetic relatedness within P. gilae was further assessed by a haplotype network 
that was generated by TCS version 1.21 using the default settings (e.g., 95% connec- 
tion limit) and fixing the connection limit at 90 steps (Clement et al. 2000). NDI 
sequence divergences (maximum composite likelihood) within and between lineages 
were calculated using MEGA6 (Tamura et al. 2013), with standard errors estimated 
by 1,000 bootstrap replications with pairwise deletion of missing data. Structuring of 
variation among lineages was evaluated by an AM OVA using Arlequin 3.5 (Excoffier 
and Lischer 2010). 

Types and other voucher material were deposited in the National Museum of 
Natural History (USNM) collection. Specimens of P. gilae from the Bell Museum 
of Natural History (BellMNH) were also examined during the course of this study. 
Series of large adults («>10) were used for shell measurements. Whorl counts refer to 
the entire shell. Sexual dimorphism in shells, which is occasionally observed in Pyr- 
gulopsis (Taylor 1987), could not be quantified owing to small sample sizes. The total 
number of shell whorls was counted (WH) for each specimen; and the height and 
width of the entire shell (SH, SW), body whorl (HEW, WBW), and aperture (AH, 
AW) were measured from camera lucida outline drawings using a digitizing pad linked 
to a personal computer (see Hershler 1989). In addition, three ratios were generated 
from the raw data (SW/SH, HBW/SH, AH/SH). Descriptive statistics were gener- 
ated using SYSTAT FOR WINDOWS 11.00.01 (SSI 2004). T-tests (two-tailed) of 
differences among shell variables were conducted using an on-line calculator (http:// 
in-silico.net/tools/statistics/ttest); data for type material of P. gilae were from Taylor 
(1987, table 11). Penial variation was described from series of adult specimens that 
were relaxed with menthol crystals and fixed in dilute formalin prior to preservation 
in 70% ethanol. Descriptive penial terminology is from Taylor (1987) and Hershler 
(1994, 1998). Variation in the number of cusps on the radular teeth («=5) was assessed 
using the method of Hershler et al. (2007). 

We used a conservative, evolutionary lineage concept in describing new species 
only for those snails that are morphologically diagnosable as well as phylogenetically 
independent and substantially divergent genetically (Hershler et al. 2007). Inasmuch as 
the principal goal of our paper was to delimit species, we provide only brief taxonomic 
descriptions which focus on those aspects of morphology that have proven most useful 
in previous such studies of Pyrgulopsis (Taylor 1987, Hershler 1994, Hershler 1998). 
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Results 

Sixteen (16) NDI haplotypes of P. gilae were detected, 1 1 of which were restricted to 
single populations (Table 2). The others were shared by pairs of populations along the 
lower East Fork (haplotype II), upper East Fork (haplotypes VII, IX, XI) and Middle 
Fork (haplotype XIII) Gila River. Six (6) samples each contained a single haplotype 
(G3, G4, G5, G7, G13, G14). The TCS analyses (not shown) recovered three well dif- 
ferentiated haplotype groups composed of specimens from along the lower East Fork 
and main stem Gila River (clade I), Middle Fork Gila River (II), and the upper East 
Fork Gila River (III). These groups differed from each other by 3.7-8.7% sequence 
divergence; variation within groups was minor (Table 3). The AM OVA indicated that 
most of the detected variation (91.7%) was partitioned among these groups; varia- 
tion within populations, and among populations within the groups was much smaller 
(1.35, 6.93%) but nonetheless was significant (Table 4). The three previously reported 
clades (I-III; Liu et al. 20 1 3) were similarly recovered in Bayesian analyses of both the 
NDI dataset, and the combined COI + NDI dataset (Fig. 2). Based on the genetic 
evidence of distinctiveness and the diagnosable shell and penial characters that are 
detailed below we recognize two of these lineages as new species which are described 
herein (clade II as P. marilynae, clade III as P. similis) and restrict P. gilae (clade I) to its 
originally circumscribed geographic range. 



Table 2. Frequency distribution of NDI haplotypes detected in P gilae. ?2=sample size. 



Haplotype (specimen code) 


Sample 




Gl 


G2 


G3 


G4 


G5 


G6 


G7 


G8 


G9 


GIO 


Gil 


G13 


G14 


I (GIA) 


2 


























II (GIC) 


1 




3 






















III (G2B) 




2 
























IV (G2D) 




1 
























V (G14B) 


























4 


VI (G4B) 








3 




















VII (G5A) 










4 


1 
















VIII (G6B) 












4 
















IX (G7A) 














8 










4 




X (G9A) 


















2 










XI (G9C) 


















2 


2 








XII (GlOA) 




















2 








XIII (G8A) 
















5 






1 






XIV (G8D) 
















1 












XV(GllF) 






















2 






XVI(GllH) 






















2 






n 


3 


3 


3 


3 


4 


5 


8 


6 


4 


4 


5 


4 


4 
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Table 3. Mean NDI sequence divergence (maximum composite likelihood) within and among B gilae 
clades. I: Gl, G2, G3, G14; II: G8, Gil: III: G4, G5, G6, G7, G9 GIO, G13. 





Glade I 


Clade II 


Clade III 


Glade I 


0.003+/-0.002 






Glade II 


0.037+/-0.011 


0.002+/-0.001 




Glade III 


0.087+/-0.021 


0.075+/-0.019 


0.006+/-0.003 



Table 4. Genetic differentiation among P gilae clades based on NDI sequences. Sub-groups=(Gl, G2, G3, 
G14), (G8, Gl 1), and (G4, G5, G6, G7, G9, GIO, G 13). Asterisked (t values are highly significant (P<0.001). 





df 


variance components 


% of variation 


•I* statistic 


Among groups 


2 


13.95 


91.72 


0.91* 


Among populations within groups 


10 


1.05 


6.93 


0.84* 


Within populations 


43 


0.21 


1.35 


0.99* 




P. gilae 
(clade I) 



P. marilynae 
(clade II) 



P. similis 
(clade III) 



mimbres" 



P. deserta 



P. inopinata 
P. plicata 
P. kolobensis 



1 00 I 

I P Unl 



■ R trivialis 

P. glandulosa 

P. sola 



-T 



P. owensensis 
P. perturbata 



Figure 2. Bayesian tree based on the combined (COI, NDI) dataset. Posterior probabilities for nodes are 
given when >95%. Specimen codes are from the Table 1. 
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Systematic descriptions 

Family Hydrobiidae 

Subfamily Nymphophilinae 

Genus Pyrgulopsis Call and Pilsbry, 1886 

Pyrgulopsis marilynae Hershler, RatclifFe, Liu, Lang and Hay, sp. n. 

http://zoobank.org/A64l736C-650D-4649-B8AD-5A012AFB3396 
Figs 3, 4A-B 

Pyrgulopsis gilae (clade II). — Liu et al. 2013. 

Types. Holotype, USNM 1 135068 (a dry shell), spring 0.48 km north, 0.48 km west 
of Jordan Canyon, Catron County, New Mexico, 33.2909°N, 108.2681°W, 1 October 
2009, Michelle Christman. Paratypes, USNM 1231474 (from same lot). 

Referred material. NEW MEXICO. Catron County: USNM 1 123432, USNM 
1123588, spring 0.8 km north, 0.64 km west of Jordan Canyon (33.2889°N, 
108.2683°W), USNM 1 135067, spring 0.97 km north, 0.64 km west of Jordan Can- 
yon (33.2924°N, 108.2696°W), USNM 883175, Jordan Hot Spring (33.2927°N, 
108.2692°W). 

Diagnosis. Distinguished from P. gilae and the species described next {P. simi- 
lis) by its narrower shell (mean shell width/shell height 0.613 vs. 0.682, t=-9.6588, 
df=36.2176, /'<0.0001, «=30 for P. gilae; 0.613 vs. 0.734, t=-l6.36l7, df=18.9656, 
P<0.0001, n=\ \ for P. similis), more pronounced whorl shoulders, and broad overlap 
of the ventral surface of the penis by the terminal gland (probably reflecting fusion 
with a distal ventral gland). Further differs from P. gilae in its smaller size (mean shell 
height2.77vs. 3.47 mm, t=-11.3848, df=21.9544, /'<0.0001) and (basal) extension of 
the outer penial gland to mid-line or left edge of penis. Further differs from P. similis in 
its larger size (mean shell height 2.77 vs. 2.36 mm, t=7.3691, df= 1 5.370 1, /'<0.0001), 
smaller number of dorsal glands on the penis, and larger size of the terminal and ven- 
tral glands on the penis. 

Description. Shell (Fig. 3A— B) narrow-conic, whorls 4.5—5.0. Teleoconch whorls 
convex, shoulders narrow, angular, sutures impressed. Aperture ovate, angled above, pa- 
rietal lip complete, usually slightly disjunct, umbilicus narrow. Outer lip thin, orthocline. 

Operculum (Fig. 3C-D) as for genus; edges of last 0.5 whorl weakly frilled on 
outer side; portion of muscle attachment margin thickened on inner side. Radula (Fig. 
3E-G) as for genus; dorsal edge of central teeth concave, lateral cusps four-five, basal 
cusp one. Lateral teeth having two-three cusps on both inner and outer sides. Inner 
marginal teeth with 14-20 cusps, outer marginal teeth with 17-22 cusps. Radula data 
are from USNM 1135067. 

Penial filament and penial lobe about equal in length (Fig. 4A-B). Filament having 
two (penial) glands on dorsal surface; inner gland shorter. Outer penial gland curving 
to mid-line (10/24 specimens) or left edge of penis (14/24 specimens), the latter con- 
dition probably represents fusion with a gland on the left edge (Dg2). Terminal gland 
elongate, horizontal, broadly overlapping ventral surface of penis. Dorsal surface of penis 
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Figure 3. Shells, opercula and radula, P. marilynae n. sp. A Holotype, USNM 1135068 B Shell, USNM 
1135067 C, D Opercula (outer, inner sides), USNM 1135067 E Portion of radular ribbon, USNM 
1 135067 F Central teeth, USNM 1 135067 G Lateral and inner marginal teeth, USNM 1 135067. Scale 
bars A, B 1.0 mm; C, D 100 |^m; E 20 |.im, F, G 10 [im. 



Table 5. Shell parameters for P. marilynae. Measurements are in mm. 





WH 


SH 


sw 


HBW 


WBW 


AH 


AW 


SW/SH 


HBW/SH 


AH/SH 


Holotype, USNM 1 135068 




4.75 


2.99 


1.78 


2.12 


1.56 


1.22 


1.15 


0.60 


0.71 


0.41 


USNM 12231474 (k=10) 


Mean 


4.60 


2.77 


1.70 


2.01 


1.51 


1.16 


1.6 


0.61 


0.73 


0.42 


S.D. 


0.13 


0.15 


0.08 


0.10 


0.07 


0.05 


0.05 


0.01 


0.02 


0.0 


Range 


4.50^.75 


2.51-3.06 


1.53-1.83 


1.83-2.15 


1.37-1.62 


1.04-1.22 


0.95-1.12 


0.60-0.64 


0.69-0.76 


0.40-0.44 



having gland along right edge of lobe (Dg3) and 2-3 additional glands (22/24 specimens); 
one specimen did not have any additional glands and one specimen had four additional 
glands. Ventral gland positioned near centrally. Penial data are from USNM 11 35067. 



78 



Robert Hershler et d. I ZooKeys 429: 69-85 (2014) 



Etymology. The specific epithet is a patronym honoring Marilyn Myers (United 
States Fish and Wildlife Service, retired) for her dedicated efforts to survey Pyrgulopsis 
habitats in the upper Gila River basin. 

Distribution. A series of seeps and springs along the north side of short reach (ca. 
0.25 km) of the Middle Fork Gila River just below Jordan Hot Spring (Fig. 1). The 
type locality is a seep wall which is the lower-most occurrence of P. marilynae along the 
Middle Fork Gila River; the water temperature at this site was 25°C on 1 October 2009. 

Remarks. Pyrgulopsis marilyane was resolved as sister to P. gilae (100% posterior 
probability) in the molecular phylogenetic analysis (Fig 2). The apparent fusion of the 
terminal and distal ventral glands of the penis that characterizes this species (in part) 
was previously reported for P. sadai (Hershler 1998, fig. 391). The sample attributed to 
Jordan Hot Spring (USNM 883175) may have been collected instead from a closely 
proximal spring as P. gilae has not been found at the former locality during recent 
surveys (USFWS 201 lb). 

Pyrgulopsis similis Hershler, RatcIifFe, Liu, Lang and Hay, sp. n. 

Figs 4C-D, 5 

Pyrgulopsis gilae. — Hurt 2004 (in part; Wall Lake population). 
Pyrgulopsis gilae (clade III). — Liu et al. 2013. 

Types. United States: Holotype, USNM 1135064 (a dry shell), spring along Beaver 
Creek, ca. 0.29 km north and 0.4 km west of confluence with Taylor Creek, Catron 
County, New Mexico, 33.3405°N, 108.1097°W, 21 May 2009, BKL and Marilyn My- 
ers. Paratypes, USNM 1135065, 1231475 (from same lot). 

Referred material. NEW MEXICO. Catron County: USNM 854684, USNM 
1135057, USNM 1123589, USNM 1135058, USNM 1135059, Fall Spring, 
1.61 km north, 0.56 km east of Burnt Corral Canyon (33.294°N, 108.1302°W), 
USNM 1123590, hillside seep 1.61 km north, 0.97 km east of Burnt Corral Canyon 
(33.2951°N, 108.1268°W), USNM 854685, USNM 1123594, USNM 1135060, 
USNM 1135061, seepage along Taylor Creek, ca. 0.32 km south, 0.93 km west 
of Wall Lake dam (33.3457°N, 108.0904°W), USNM 854683, USNM 1123592, 
USNM 1135062, USNM 1135063, spring along Taylor Creek, ca. 0.81 km north, 
1.13 km east of Wall Lake Dam (33.358rN, 108.0673°W), USNM 854682, USNM 
1123593, NM: Catron Co., seepage along Taylor Creek, ca. 50 m west of Whitetail 
Canyon (33.3613°N, 108.0576°W). 

Diagnosis. Difi^ers from P. gilae in its smaller size (mean shell height 2.36 vs. 3.47 
mm, t=-22.7297, df=36.4071, /'<0.0001, n=30 for P. gilae), larger number of glands 
on the dorsal surface of the penis, frequent extension of outer penial gland and/or Dg2 
to the mid-line of the penis, and smaller size of the terminal and ventral glands on the 
penis. Contrasted with P. similis above. 
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Figure 4. Penes (dorsal, ventral surfaces). A, B P. marilynae n. sp., USNM 883175 C, D P. similis n. 
sp., USNM 1135065 E, F P. gilae, BellMNH 20898. Scale bars A-F 200 (jm. Dg2 dorsal gland along 
left edge Dg3 dorsal gland along right distal edge Ipg inner (left) penial gland Opg outer (right) penial 
gland Tg terminal gland Vg ventral gland. 



Description. Shell (Fig. 5A-B) ovate- to narrow conic, whorls 3.75-4.50. Tel- 
eoconch whorls medium convex, narrowly shouldered. Aperture pyriform, parietal 
lip complete, usually adnate, sometimes slightly disjunct, umbilicus small. Outer lip 
thin, orthocline. 

Operculum (Fig. 5C-D) as for genus; edges of last 0.5 whorl frilled on outer 
side; inner side near smooth. Radula (Fig. 5E-G) as for genus; dorsal edge of central 
teeth concave, lateral cusps four-six, basal cusp one. Lateral teeth having two-three 
cusps on inner sides and two-four cusps on outer sides. Inner marginal teeth with 
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Figure 5. Shells, opercula and radula, P. similis n. sp. A Holotype, USNM 1 135064 B Shell, USNM 
854684 C, D Opercula (outer, inner sides), USNM 1135065 E Portion of radular ribbon, USNM 
1 135065 F Central teeth, USNM 1 135065 G Lateral and inner marginal teeth, USNM 1 135065. Scale 
bars A, B 1.0 mm; C, D 200 ^m; E 20 |im, F, G 10 fim. 



15-20 cusps, outer marginal teeth with 16-25 cusps. Radula data are from USNM 
1135059, USNM 1135064. 

Penial filament longer than lobe (Fig. 4C— D). Filament having two (penial) glands 
on dorsal surface; inner gland shorter. Outer penial gland sometimes extending (basal- 
ly) to mid-line (4/30 specimens) or left edge (7/30 specimens); Dg2 sometimes curv- 
ing (basally) to mid-line (11/30 specimens). Terminal gland transverse, rather small. 
Dorsal surface of penis having gland along right edge of lobe (Dg3) and 3-7 additional 
glands (30/30 specimens) which form long, slightly oblique strips. Ventral gland small, 
positioned near centrally; second gland rarely present (4/30 specimens). Penial data are 
from USNM 1135065. 



Taxonomic revision of the Pyrgulopsis gilae (Caenogastropoda, Hydrobiidae)... 81 



Table 6. Shell parameters for P similis. Measurements are in mm. 





WH 


SH 


sw 


HBW 


WBW 


AH 


AW 


SW/SH 


HBW/SH 


AH/SH 


Holotype, USNM 1135064 


4.25 


2.52 


1.73 


1.94 


1.48 


1.24 


1.10 


0.69 


0.77 


0.49 


USNM 1231475 


n=n) 


Mean 


4.05 


2.36 


1.73 


1.89 


1.50 


1.22 


1.04 


0.73 


0.80 


0.52 


S.D. 


0.10 


0.10 


0.09 


0.07 


0.07 


0.06 


0.04 


0.02 


0.02 


0.02 


Range 




2.16-2.50 


1.62-1.89 


1.78-2.01 


1.41-1.61 


1.14-1.30 


0.98-1.09 


0.71-0.76 


0.78-0.83 


0.49-0.55 



Distribution. Springs along a short reach (ca. 1 0 km) of the East Fork Gila River 
from just above Wall Lake to slightly above the mouth of Burnt Corral Canyon (Fig. 1). 
The type locality is a spring brook (ca. one m wide and 0.25 m deep) that discharges at 
the base of the canyon wall along the east side of Beaver Creek; the water temperature 
at this locality was 22.1°C on 21 May 2009. The flow at this locality is augmented by 
numerous small seeps. 

Etymology. The specific epithet is an adjective referring to the close resemblance 
between this species and both P. gilae and P. marilynae. 

Remarks. Pyrgulopsis similis was resolved as sister to the clade composed of P. 
marilynae and P. gilae (100% posterior probability) in the Bayesian analysis of mo- 
lecular data (Fig. 2). 

Pyrgulopsis gilae (Taylor, 1987) 

Fig 4E-F 

Fontelicella gilaeTaylor 1987: 16-18, fig. 7,tables 11-13 (springs on north side of East Fork 
of Gila River, center of sec. 3, Tl 3S, Rl 3W, unsurveyed. Grant County, New Mexico). 
Pyrgulopsis gilae. — Hershler 1994: 36-38, figs 15a-c, 46c (new combination). 
Pyrgulopsis gilae. — Hurt 2004 (in part; Gila I-III populations). 
Pyrgulopsis gilae (clade I). — Liu et al. 2013. 

Types. Holotype, LACM 2214; paratypes, BellMNH 20898, BellMNH uncat., UTEP 
10054, USNM 854087 (from same lot as holotype). 

Referred material. NEW MEXICO. Grant County: USNM 1135050, USNM 
1135052, spring ca. 1.29 km mile north, 0.55 km west of confluence of East Fork 
Gila River and Black Canyon (33.1864°N, 108.1675°W), USNSM 1 123426, USNM 
1135055, USNM 1135056, spring ca. 1.53 km north, 2.38 km east of State Route 
527 bridge crossing (33.1946°N, 108.1804°W). 

Other material examined. NEW MEXICO. Grant County: topotypes, USNM 
1004620, USNM 1135043, USNM 1135044, spring ca. 1.53 km north, 2.90 km 
east of State Route 527 bridge crossing (33.1917°N, 108.1742°W), BellMNH uncat., 
USNM 873211, USNM 1068942, "Alum Hot Spring," ca. 1.93 km south, 0.16 km 
west of State Route 527 bridge crossing (33.1618°N, 108.208rW). 
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Figure 6. Shells, P. gilae. A USNM 1135052 B USNM 854574. Scale bars A, B, 1.0 mm. 



Distribution. Several groups of springs in the lower reach of the East Fork Gila 
River (below the mouth of Black Canyon) and a single spring along the Gila River ca. 
2 km below the East Fork confluence (Fig. 1). 

Remarks. Examination of the large series of penes that Taylor scored for this spe- 
cies (BellMNH 20898, BellMNH uncat.) indicated that neither the outer penial gland 
nor Dg2 extends appreciably onto the dorsal surface of the penis (Fig. 4E-F; also see 
Taylor 1987, fig. 7b-c) in contrast with P. marilynae and P. similis. Specimens from 
the two new populations (Fig. 6A), which are closely proximal to the type locality, and 
the disjunct "Alum Spring" population (Fig. 6B) conformed to P. gilae in all morpho- 
logical details. Pyrgulopsis gilae co-occurs sympatrically with P. thermalis (Taylor) at 
several localities (Taylor 1987). 



Discussion 

The results of this study provide additional evidence that the current taxonomy of Pyr- 
gulopsis in some cases masks cryptic species diversity. Our previous revision of widely 
ranging P. micrococcus revealed this taxon to be a polyphyletic composite of five species, 
three of which were undescribed (Hershler et al. 2013). Here we have shown that, in 
contrast, P. gilae (in the broad sense) is a monophyletic species complex diagnosed by 
a unique penial character — the presence of two glandular strips on the dorsal surface 
penial filament. (Note that P. merriami [Pilsbry and Beecher] , which is distributed in 
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an isolated basin in southeastern Nevada, has a somewhat different pattern consist- 
ing of two glands on the dorsal and one gland on the ventral surface of the filament; 
Hershler 1994). These disparate findings underscore the complexity of and taxonomic 
challenges posed by the Pyrgulopsis radiation, which is characterized by endemism on 
very fine geographic scales and extensive morphological homoplasy (Hershler 1994, 
Liu and Hershler 2005). 

The delineation of cryptic species complexes often has important consequences 
for conservation (Bickford et al. 2007). Pyrgulopsis gilae was recently removed from 
the federal candidate list for listing as endangered or threatened in part owing to the 
discovery of populations from along the East Fork and Middle Fork of the Gila River 
(USFWS 201 la, USFWS 20 lib) that are assigned herein to two new species. The re- 
sulting restriction of P. gilae to its originally circumscribed geographic range — several 
groups of springs in the lower reach of the East Fork Gila River (below the mouth of 
Black Canyon) and a single spring along the Gila River ca. 2 km below the East Fork 
confluence — suggests the conservation status of this species should be re-visited by the 
USFWS. The narrow endemism of the two new species suggests that these may also 
merit consideration for possible listing by the USFWS. 

Our findings also underscore the need for additional field surveys to further deline- 
ate the occurrences of Pyrgulopsis in New Mexico and to supplement the recent mono- 
graph by Taylor (1987). Large portions of the Gila River and other drainage basins in 
the state have yet to be carefully searched for these tiny animals. 
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